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Abstract 

I develop a diagrammatic method for calculating chiral logarithms in the 
quenched approximation. While not rigorous, the method is based on physically 
reasonable assumptions, which can be tested by numerical simulations. The main 
results are that, at leading order in the chiral expansion, (a) there are no chiral 
logarithms in quenched f n , for m u = m^; (b) the chiral logarithms in Bk and 
related kaon B-parameters are, for = m s , the same in the quenched approxi- 
mation as in the full theory; (c) for and the condensate, there are extra chiral 
logarithms due to loops containing the rf, which lead to a peculiar non- analytic 
dependence of these quantities on the bare quark mass. Following the work of 
Gasser and Leutwyler, I discuss how there is a predictable finite volume depen- 
dence associated with each chiral logarithm. I compare the resulting predictions 
with numerical results: for most quantities the expected volume dependence is 
smaller than the errors, but for By and Ba there is an observed dependence 
which is consistent with the predictions. 
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1 INTRODUCTION 



The approximate chiral symmetry of QCD provides strong constraints on the form of am- 
plitudes involving pseudo-Goldstone bosons. A classic example is that the pion scattering 
amplitudes at threshold are related to the pion decay constant |l| 

2 

Tr(s = 4ml,t = u = 0) = c I -f + O(ml) . (1) 

J n 

The cj are known numerical constants which depend on the isospin / of the pion pair. The 
corrections of 0(m*) fall into two classes:^ those that are analytic in m 2 , coming from 
non- leading terms in the chiral Lagrangian ||; and the non-analytic "chiral logarithms" [| 
proportional to m^.ln{m 7r ). These logarithms come from infrared divergences in pion loops 
[||. [I| . They have the important property that they are proportional to the lowest order term 
in the chiral expansion. For example, in the limit m n = m^ = 0, f n has the expansion ||. [2|] 

U = f{l-\L(m K ) + cm 2 K ) . (2) 

where the chiral logarithm is 

(I use the normalization such that f n = 93 MeV.) Chiral symmetry alone does not determine 
/, A or c (in fact c can be absorbed into A), but does fix the coefficient of the chiral log L. 

In the continuum, chiral logs are useful as a guide to the size of the corrections to the 
leading order chiral behavior. If the logs are large when one takes a reasonable scale A ~ m p , 
then there is reason to worry about the reliability of leading order chiral predictions. This 
approach has been taken for example in Refs. Q and ||. 

On the lattice, chiral logs are potentially of greater importance, because they allow one 
to test whether the contribution of pion loops is being correctly described. These loops give 
rise not only to chiral logs, but also to the pion cloud surrounding hadrons. This cloud is 
an important part of the structure of hadrons, affecting charge radii and other form factors. 
The obvious way to search for the chiral log, by studying the dependence of amplitudes on 
m 2 or m 2 K , is not very useful, for it is difficult to disentangle a logarithm from a power series. 
It is better to study the volume dependence of the amplitude. As explained by Gasser and 
Leutwyler |]], associated with each chiral log is a known volume dependent correction. For 
the particular case of the masses, the correction can also be obtained from a general formula 
given by Liischer 0. 

Chiral logarithmic corrections have been calculated for a number of a quantities in QCD. 
Those that are of interest here are the pseudo-Goldstone boson masses and decay constants, 
the condensate, and Bk- I have extended these calculations to By and Ba- These results 
can be converted into predictions of finite volume dependence for lattice simulations of full 

2 Here and in the following I often use "pion" to generically to refer to all of the pseudo-Goldstone bosons, 
i.e. the pions, kaons and eta's. In the present example, 0(jn%) includes the possibility of Oijn^m 2 ^). 

3 For baryon properties some of the non-analytic terms are not logarithms, but are proportional to (m 2 ) 3 / 2 
0. I do not consider baryon properties here. 
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QCD, and I give some examples below. At present, however, most lattice calculations of 
these quantities use the so-called "quenched" approximation, in which internal quark loops 
are not included. The lack of internal quark loops means that many, though not all, of the 
pion loops present in full QCD are absent in the quenched approximation. The main focus 
of this article is to understand how to adapt the calculation of chiral logarithms in full QCD 
to the quenched approximation. This is done here for the quantities just mentioned, but the 
method can be extended to many other quantities. These results allow one to predict the 
finite volume dependence of quenched matrix elements. 

In addition to their diagnostic power, chiral logs provide an estimate of the size of the 
error introduced by the quenched approximation. In most quantities the chiral logs have a 
different coefficient in quenched and full theories. For example, it is shown below that there 
are no leading order chiral logs in /„. in the quenched approximation,^ in contrast to the result 
for full QCD, Eq. (0). Taking A = m p , the chiral log in full QCD gives a 6% contribution 
to f n . This suggests that the absence of chiral logs in the quenched approximation should 
have an effect on f n of this size. This is a suggestion, and not a calculation, because the 
quantities / and c need not be the same in full and quenched QCD. Nevertheless, even a 
rough guide is useful since we have so few tools for studying the effect of quenching. 

The idea of calculating chiral logarithms in the quenched approximation was first sug- 
gested by Morel ||, who was inspired by numerical evidence of non-analytic terms in m^/m q 
(Ref. ||). In order to study the quenched approximation, he introduced scalar-quark ghost 
fields which, upon functional integration, produce a determinant which cancels that from 
the integration over quarks. He showed how to calculate chiral logs in the combined strong 
coupling, large dimension limit, both for color gauge groups SU{2) and £77(3). In Ref. [1(J 
I corrected errors in Morel's result, finding that, for three colors, there were no quenched 
chiral logs in f n or m n in the large d, strong coupling limit |] More importantly, I showed 
how the presence or absence of chiral logs could be understood by studying quark diagrams. 
It is the extension of this "diagrammatic method" that I use in this paper. 

As will become clear in the following, this method is neither rigorous nor systematic. The 
lack of rigor is common to all calculations of quenched chiral logs. Such calculations use chiral 
perturbation theory, the validity of which is predicated on the existence of a well defined 
underlying field theory. Quenched QCD is not a well defined field theory. In particular, 
quenched amplitudes are neither unitary nor analytic. Nevertheless, as I discuss in Sections 
|2| and ^, it is reasonable to assume both that chiral logs are present, and that they can be 
calculated using an adaptation of chiral perturbation theory. 

The diagrammatic method is also not systematic. It amounts to a set of rules for dis- 
carding some of the chiral perturbation theory diagrams, and for changing the masses of 
the pions that appear in some of the loops of the diagrams that are kept. These rules are 

4 Strictly speaking this is only true if m u — ma, as discussed below. 

5 The reader may be confused by two points upon consulting Ref. JhJ. (1) A class of diagrams involving 
rf loops, which are discussed at length in this article, were overlooked in Ref. ||l0| . This does not affect main 
calculation of that reference, because these diagrams are not present in the large d, strong coupling limit 
||. It does, however, mean that the conclusion drawn in that paper, namely that there are no leading order 
chiral logs in quenched f n and in general, is false. This error is rectified here. (2) There are chiral logs 
in /„• and m n for color SU(2), even in the large d, strong coupling limit. These may explain the numerical 
results of Ref. ||. SU(2) is special, however, since there are Goldstone baryons as well as Goldstone pions. 
In this paper I consider only three or more colors. 
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developed in an ad hoc way, quantity by quantity. I do not know how to write them down 
in general. Bernard and Golterman [II] have recently outlined a more systematic approach 



in which chiral perturbation theory is extended to include states containing Morel's scalar 
ghosts. Their method appears to give the same answers for the quantities that I calculate 
here. I prefer to use the diagrammatic method, however, because it displays the underlying 
physics very clearly. For more complicated calculations Bernard and Golterman's method is 
likely to be superior. 

A difficulty with calculations of quenched chiral logarithms stems from the fact that one 
must treat the flavor singlet rj as a pseudo-Goldstone boson. This is in contrast to full QCD, 
where the rf is not a pseudo-Goldstone boson because of the axial anomaly. It turns out the 
contributions of rf loops do not satisfy the usual power counting rules of chiral perturbation 
theory, in which each new loop brings with it a power of m 2 . Thus corrections proportional to 
powers of ln(m,r) are obtained from diagrams with any number of loops. This is a potentially 
disastrous, for it might undermine the application of quenched chiral perturbation theory 
to all quantities. As a first attempt to understand the issue, I have summed the leading 
logarithms, and find that they can be absorbed into a change of definition in the quark 
mass. The only observable effects are that m n , (if)il>), and related quantities should have a 
weak finite volume dependence. More study of this issue is required. 

Fortunately, it turns out that there are a class of "safe" quantities which are insensitive, 
at leading order, to rj' loops. Examples are the decay constants and kaon B— parameters for 
degenerate quarks. I concentrate on these quantities for most of this paper. 

Bernard and Golterman have suggested a more ambitious approach They point 

out that, although rf loops are not parametrically of higher order in chiral perturbation 
theory, they are non-leading in a combined chiral and large N c expansion. Furthermore, 
their contribution is numerically small. Thus, from a practical point of view, it may make 
sense to do a loop expansion, and they have begun such a program. 

The work presented here was inspired by numerical results showing volume dependence in 
By and Ba jl^l - It turns out that these results are described reasonably well by predictions 



given here. Present calculations of other quantities are not yet accurate enough to test the 
predictions. 

The outline of this paper is as follows. In the next Section I recall some useful facts 
about the chiral Lagrangian, and discuss its extension to the quenched approximation. In 
Section [| I introduce the method of calculating chiral logs in the quenched approximation 
using the example of decay constants. Section |] extends this method to Bx and the related 
parameters By and B^. In Section [5] I discuss the contribution of rf loops, and in Section 
|^ I explain how to predict the leading finite volume dependence from the results for chiral 
logs. These predictions are confronted with numerical data in Section [?]. Section § contains 
some conclusions. 

Brief discussions of some of the results of this paper have been given previously in Refs. 



2 CHIRAL LAGRANGIAN 

Chiral logarithms in full QCD are most simply calculated using the chiral Lagrangian. In 
this section I first collect some standard results, and then explain my assumptions for how 
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these are modified in the quenched approximation. I use the notation of Ref . |16| , and work 
in Euclidean space. The Lagrangian is 



C = \f Tr (d v Zd v tf) - \fp, Tr (M(S + E*)) + 0(p\ Mp 2 , M 2 ) , (4) 

where £ = exp(2iii a T a / f), the n a being the pseudo-Goldstone bosons fields. The group 
generators T a are normalized such that Tr(T a Ti,) = ^5 a b. M is the quark mass matrix, in 
terms of which the pion mass matrix is (m 2 ) a b = 4/i Tr(MT a Tt>). The form of £ is identical 
for any number of flavors. Flavor dependence enters only through the number of pions (N 2 — 1 
for flavor SU(N)), and through the constants / and /i. 

In QCD the mass matrix is diag(m M , m<j, m s ). I will assume that m u = which is a 
good approximation for the amplitudes I discuss, as the terms that I drop are suppressed 
by {rn u — m^)/m s compared to those that I keep. In this isospin symmetric limit the flavor 
off-diagonal eigenstates are the 7r ± , with m 2 = and the kaons K ± , K , Kq with m 2 K = 

fj,(m s + m d ), while the flavor diagonal eigenstates are the n° (T a = T 3 ), which is degenerate 
with the 7r ± , and the 77 (T a = T 8 ), with m 2 = |/i(4m s + 2m d ). 

The left- and right-handed quark currents can be written in terms of £ fields 

L, a = qT al ^(l+ l5 )q — \f Tr(T a W^) , (5) 
R,a = qTal,U 1 -^)q — i/ 2 Tr(T a S^S) . (6) 

The Euclidean space Dirac matrices are hermitian and satisfy {7^, 7^} = 5^. The expansion 
of Lfj, in terms of the pion field is 

Lya = -if (Tr(T a «V) - j Tr(T a [«V, A) ~ Tr(T a [7r[7r, d^}}) + 0(tt 4 )^ . (7) 

has the same expansion except that it — > —71. The pion decay constant is defined by 

(0 1 (7^7^759 K) = (0\L^a-Rf,a\n) = Ph/tv ■ (8) 

Comparing with Eq. ([?]) we see that f n = fa,t lowest order in chiral perturbation theory. 

The chiral Lagrangian can be extended to include the flavor singlet pseudoscalar meson 
[0]. Throughout this paper I refer to this meson as the rj', despite the fact that the physical 
rf is not exactly a flavor singlet. I use the label "pion" to refer to the non-singlet pseudo- 
Goldstone bosons, not including the rf. To include the rf in the chiral Lagrangian one begins 
by defining a U(N) matrix 

U = £e^° ; 0o = 2 V '/(V2Nf) , (9) 

in terms of which the extended Lagrangian can be written 

C = +i/V 1 (?/)Tr[(^f/,%f/ t ]-^^^ 

+Vo(v') + V 5 (r ] WdvV' ■ (10) 

(The definitions of the potentials differs slightly from those of Ref. 0.) The potential are 
arbitrary aside from the following conditions: V , V\, and V 5 are real, and even functions of 



4 



their arguments; V2 satisfies V^t/)* = V^— 7/); and Vi(0) = V^O) = 1. The arbitrariness 
implies that there are very few constraints on the couplings of the rj' '. 

This Lagrangian simplifies when N c , the number of colors, is taken to infinity. This 
limit is of interest here because it is similar to the quenched approximation. From the 
analysis of Ref. 0, it can be seen that, when N c — > 00, one can set V± — V2 — 1, and 
Vo = V5 = 0. (Vq does contribute a constant term proportional to A^, but this only affects 
the vacuum energy.) In other words, the Lagrangian has the same form as the original chiral 
Lagrangian, Eq. @, except that E is replaced by U. Dependence on N c enters only through 
the decay constant, which is proportional to ^fW c . The symmetry group has thus enlarged 
to U(N)l x U(N)r, with the rf being the additional pseudo-Goldstone boson. The physical 
reason for this simplification is that the terms in Eq. ( |T0D which differentiate the rf from 
the pions involve intermediate gluons, and are suppressed by powers of 1/N C . 

Expanding the large N c Lagrangian in terms of E and rj' one finds 

£(iV c ^oo) = £ + £ m , (11) 

£0 = \f (rt [dj&jf] + ^i/0 M t/J , (12) 

Cm = -|/V(cos(^ o )'L:[M(E+E t )]+iBin(0o) , &[M(E-E t )]) . (13) 
Three features of this result are important in the following. 

• The rf couples to pions through the mass term C m , but does not couple through the 
derivative term Co- This is not true at finite N c because the potential V\ will contain 
even powers of the rf field. 

• The masses of the pseudo-Goldstone pions differ in general from those at N c = 3. 
Diagonalizing the mass matrix, one finds that the rj f , rj and tt° mix. The eigenstates 
have flavor compositions ss, dd and uu, with masses 2/im s , and 2fim u respec- 
tively. For rrid = m u one can use the usual ir° and rj linear combinations, since they 
are degenerate. Only if all three quarks are degenerate, however, is the flavor singlet 
rj' an eigenstate. 

• The non-singlet chiral currents that follow from the Lagrangian are the same as those 
of Eqs. (H) and (|j), except that E is replaced by U. It is simple to see, however, that 
the rj' drops out of the expressions, so that the currents are in fact unchanged. This 
can be seen explicitly in the expansion of L^ a , Eq. (0). The first term vanishes for the 
rj' because of the flavor trace, and all subsequent terms involve commutators, which 
are zero for the rj'. 

The fundamental assumption of this paper is that the "non-hairpin" interactions involv- 
ing pions and rj's are described by a Lagrangian which has, at leading order in the chiral 
expansion, the same form as the large N c Lagrangian, Eq. (|rT|). The constants / and ii will, 
however, differ from those in both full QCD and the iV c — > 00 limit. There are additional 
"hairpin" interactions involving the rj', and possibly pions, but not involving pions alone. 
Hairpin interactions are those in which the quark and antiquark fields in at least one of 
the external rj' mesons are contracted together. Examples appear in Figs |3|(d)-(i). Before 
discussing these additional interactions, I explain the justification for the assumption just 
announced. This assumption has three parts. 
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1. Interactions involving different numbers of non-singlet pions, but no flavor singlet 
mesons, are related to one another in the same way as in full QCD. For example, 
the pion scattering lengths in the chiral limit can be expressed in terms of f n accord- 
ing to the formulae of Weinberg ffl. In full QCD the form of the chiral Lagrangian 



follows from three inputs [TT\: (a) Ward identities between correlation functions; (b) 
pion domination of low momentum amplitudes; and (c) analyticity and unitarity of the 
amplitudes. In the quenched approximation one has only the first two ingredients: the 
Ward Identities in the quenched approximation are the same as in full theories, at least 
in lattice regularized theories [Tj^ and pion dominance is a good approximation in 



practice. But the amplitudes are neither analytic nor unitary. In fact, it is only ana- 
lyticity that is important for the 0(p 2 ) terms in the Lagrangian; the amplitudes from 
these terms are not themselves unitary. What is required is that multi-pion amplitudes 
have a polynomial dependence on the kinematic invariants. This is the real assump- 
tion being made concerning the quenched approximation. For staggered fermions, it is 



supported by numerical evidence for amplitudes involving up to four pions |20j. This 
is actually sufficient for the chiral logs discussed here, which require only the two and 
four point vertices from the 0(p 2 ) terms in the Lagrangian. 

2. Non-hairpin interactions of rf mesons with pions come entirely from C m , and are re- 
lated to interactions of pions with each other by flavor factors. In the large N c limit 
this result follows because the quark diagrams are the same with or without rf mesons. 
This equivalence of diagrams is also true in the quenched approximation. In particular, 
there are no internal quark loops in either the large N c limit or the quenched approx- 
imation. Interactions involving intermediate gluons, which vanish when N c — > 00, 
involve hairpins, and are discussed below. This distinction is useful because it turns 
out that, for some quantities, such hairpin interactions do not contribute at leading or- 
der. This will become clearer in the following sections. Thus I assume that the form of 
the effective Lagrangian describing the pion-?/ interactions is the same in the quenched 
and N e — > 00 theories. The coefficients in the Lagrangians will differ, however, because 
the gluon interactions are different. 

3. The interaction of the rf with itself, through both the the kinetic and mass terms, is 
related to the self-interactions of pions by flavor factors. This applies only to non- 
hairpin interactions resulting from diagrams in which the rf does not annihilate into 
gluons. The justification for this assumption is as for the previous one: the result is 
true in the large N c limit, and the quenched diagrams involving rf mesons are related to 
those involving pions by the same flavor factors as in the large N c limit. In particular, 
the mass eigenstates are the same as in that limit. 

An important consequence of the assumption that the quenched Lagrangian has the same 
form as that for C(N C — > 00) is that the non-singlet chiral currents are also the same. As 
explained above, this means that the rf does not couple to these currents. 

Although the form of the quenched Lagrangian is similar to that for full QCD, the calcula- 
tion of loop corrections using chiral perturbation theory involves additional rules. According 
to these rules certain loop diagrams are kept, while others are discarded. The chiral pertur- 
bation theory diagrams which are kept are those corresponding to quark diagrams which do 



6 



not contain internal quark loops. The method by which one determines the correspondence 
between pion and quark diagrams is best explained by example, and is the subject of the 
following Section. 

I now return to the additional hairpin interactions involving the rf . It is here that the 
quenched and large N c theories differ: such diagrams exist for the quenched approximation, 
but vanish in the N c — > oo limit. For example, there is a "disconnected" contribution to the 
rf two-point function in which the quarks annihilate into gluons. It is this diagram, and its 
iterates including any number of quark bubbles, that give the rf its additional mass in full 
QCD. More generally, such diagrams give rise to the rf dependence in the potentials V , V\, 
V 2 and V 5 , and thus convert the Lagrangian £ + C m into one with the form of Eq. QlQf) . 

In quenched calculations one must include the contribution from rf loops since the rj' is 
light. This means that the additional vertices coming from the V{ may be needed. Thus the 
coefficients of chiral logs due to rf loops are not, in general, predicted in terms of the leading 
order coefficients.^ They depend on new parameters which even in full QCD are either poorly 
constrained by phenomenology, or not known at all. It turns out, however, that there is a 
class of quantities for which the rf does not enter at one-loop order, so that the additional rf 
interactions are not needed. This class includes f n and the kaon B-parameters, in the limit 
that all quarks are degenerate. These "safe" quantities are the main focus of this article. 

For the remaining quantities one must include rj' loops. As I explain in Section |5], these 
loops do not obey usual power counting rules of chiral perturbation theory, according to which 
each loop brings an additional power of m 2 . Instead, the leading order chiral logarithms are 
not suppressed by m 2 , and come from diagrams with any number of rf loops. These diagrams 
are of a simple form, however, and require only the additional rj' 2 inter act ions. [] These can 
be parameterized as 

= ^dWdW + ^fV , (14) 

where the first term comes V5, the second from Vq. The most important feature of is that 
ttiq is a fixed scale, of O(Aqcd), and does not vanish in the chiral limit. The parameterization 
of Eq. (14) is chosen so that, if this vertex could be iterated as in full QCD, the rj' would 

get a mass 

2 2 2/im 

<> = ™o + -^r ■ ( 15 ) 



In the approximation that A = 1, the experimental rj and rj' masses imply that m ~ 0.9 
GeV. I use these values when making estimates of the quantitative importance of rj' loops. 



6 In full QCD this problem does not arise because the rj' is heavy, and does not give rise to chiral logs. 
The uncertainties only enter through the 0(p A ) terms in the chiral Lagrangian. 

7 Strictly speaking there are also two terms involving V% which contribute. (1) A mass term proportional 
to V^"(0)f/ 2 Tr(M). Using arguments similar to those of the following Section, one finds that the Tr(M) 
corresponds to an additional quark loop. Since this loop is absent in the quenched approximation, this term 
should not be included. (2) A vertex proportional to (0)77' Ti(Mt]'tt 2 ). This vertex actually appears in 
Fig. ^(f), and does contribute to mass renormalization. By redefining the U field, however, one can absorb 
V 2 ' (0) into the constant A Thus there is no loss of generality if one excludes this term. 
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(a) (b) (c) 

Figure 1: Diagrams contributing to /„■: (a) leading order; (b) wavefunction renormalization; (c) vertex 
renormalization. The lines represent pseudo-Goldstone bosons, the box is the axial current, and the circle is 
the vertex coming from C 

3 CHIRAL LOGARITHMS IN f n AND f K 

I begin by summarizing the calculation of the chiral logarithmic corrections to the decay 
constant in a theory with N light flavors, with quark loops included. In chiral perturbation 
theory there are two one-loop diagrams, and they are shown, along with the leading order 
diagram, in Fig. 1. The circle represents the four-pion vertex obtained by expanding Eq. 
(f|). The square represents the axial current ^(L^ a —R^ a ), the expansion of which is as in Eq. 
(0), except that only terms with an odd number of pions appear. The crucial point about 
these vertices is that they are completely determined in terms of the lowest order parameters 
/ and fi. 

As explained in the previous section, the leading order result is f{^d) — f, for all flavors 
d. The corrections from Figs. 0(b) and |l](c) are multiplicative 

f(rr d ) = /(l + X(Fig.l(b)) + X(Fig.|(c)) + 0(mD) . (16) 

The non-logarithmic corrections proportional to m 2 come from 0(p 4 ) terms in the chiral 
Lagrangian. They do not give rise to finite volume dependence (see Section []), and I will 
drop them from subsequent equations. 

Figure [TJ(b) gives rise to both mass and wavefunction renormalization, but only the 
latter effects f^. To contribute to wavefunction renormalization, the vertex in Fig. |l](b) 
must contain derivatives, and thus must come from £ Q . This is an important observation to 
which I return below. The contribution of Fig. |l](b) is 

X(Fig.l(b)) = -^^Tr([T c ,[T c ,T,]]T d ) / G{k,m c ) , (17) 

6J c Jk 

where the measure is J k — j d A k / '(27r) 4 , and the integrand is the pion propagator G(k,m) = 
l/(/c 2 + m 2 ). The sum runs over the N 2 — 1 pions, and is to be done in the mass eigenstate 
basis, in which the mass of the pion of flavor T c is labeled m c . 

Figure |l|(c) renormalizes the vertex; its contribution turns out to be proportional to 
wavefunction renormalization 

X(Fig.l(c)) = -4X(Fig.|(b)) . (18) 
Thus the fractional correction to the decay constant is 

f{7rd) ~ f = -£Tr([T c , [TciTd]]Td) j l(mc) . (19) 

J c 
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Here I have defined one of the standard integrals that appear at one-loop 



h(m) = -J k G(k,m) . 



(20) 



The integral must be regulated, and I use a sharp cut-off, \k\ < A.0 The result is 



A 2 (m 2 \ 



(21) 



where L(m) is chiral logarithm defined in Eq. (|3|). The quadratic divergence can be absorbed 
by a redefinition of /. 

For QCD, the general formula gives 



in agreement with Ref. ||. 

I now turn to the quenched theory. The result of the following discussion is extremely 
simple: except for rf loops, which only contribute if the quarks are not degenerate, there are 
no leading order chiral logs in decay constants. I explain this conclusion in gory detail, as 
the insight obtained can be applied directly to more complicated calculations. 

The chiral logarithm comes from the infrared part of the loop integral, \p\ < m. The 
coefficient of the logarithm is determined by the residue of the pole at p 2 = —m 2 . In general 
this residue is itself determined by a physical, on-shell amplitude: for wavefunction renor- 
malization it is the part of the four-pion scattering amplitude proportional to the kinematic 
variable t, evaluated at threshold; for vertex renormalization it is part of the amplitude for 
the axial current to convert one pion to two pions. In other words, to calculate the coefficient 
of the logarithm one can factorize the loop in Figs. |l|(b) and (c) into two pieces: an on-shell 
pion propagator, and a physical scattering or conversion amplitude. Since these amplitudes 
are evaluated for small p, they are well described by lowest order chiral perturbation theory, 
at least for small enough m. This is why the coefficient of the log can be calculated in terms 
of only the coefficients of the lowest order Lagrangian. 

It is because of this factorization that one can hope to calculate the coefficient of the 
logarithm in the quenched approximation. As discussed in the previous section, it is reason- 
able to assume that the leading order expressions for both the scattering amplitudes and the 
chiral currents are the same in the quenched approximation as in the full theory, except for 
differences in the coefficients / and fi. Thus the coefficient of the chiral logs coming from 
each diagram is the same function of these coefficients. In other words, the calculation of 
individual chiral perturbation theory diagrams proceeds by the same method in both theo- 
ries. The only differences in the quenched calculation are that the rf can propagate in the 
loops, and those diagrams which involve internal quark loops should be discarded. 

8 As is well known, such a momentum cut-off violates chiral symmetry. This can be corrected by adding a 
quartically divergent counterterm, but this does not contribute to any of the processes considered here |2l| . 



fx 

fn 




(22) 
(23) 



in agreement with the Refs 



H, 0, ^2|. For N degenerate quarks the result is 




(24) 
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Figure 2: Quark diagram corresponding to Fig. 1(a). The lines represent quark propagators, as described 
in the text. The box is the axial current. The apex where the propagators meet is the operator which is 
used to create the pion. 

To determine what these differences mean in practice it is necessary to find the quark 
diagrams which correspond to the diagrams of chiral perturbation theory. Figure |2] shows 
the quark diagram corresponding to Fig. [3](a). The quark lines are to be interpreted as 
follows. If one were doing a lattice (or any other non-perturbative) calculation of f n one 
would calculate the diagram of Fig. ||, with the lines corresponding to quark propagators 
in a background gauge field. The result would then be averaged over the all background 
fields, using either the full measure including quark loops, or the quenched measure. The 
resulting value of f n would be correct to all orders in m^, and in particular would include 
the chiral logarithms. In fact, the lines in Fig. |2| are almost the full propagators, except 
that the propagation is restricted so that there is no long distance pion contribution. If 
the full measure is used, this restriction applies also to the quarks loops contained in the 
background field. With this restriction, the diagram does not produce chiral logarithms. It 
does contribute terms of all orders in m^, but we are interested only in the lowest order 
contribution. 

This provides a definition of the meaning of the lines in quark diagrams that is both 
vague and difficult to implement. These shortcomings are not important, however, because 
the diagrams serve only as a method for tracing the flow of the flavors of the quark and 
antiquark. This aspect of the diagrams is completely well defined. 

The quark diagrams contributing to wavefunction renormalization (i.e Fig. [l|(b)) are 
shown in Fig. [3| These diagrams are obtained as follows. First, draw all possible quark 
diagrams for the pion scattering amplitude in which at least two of the pions have the flavor 
under consideration. One of these two is the external pion, the other couples to the axial 
current. Next, join the remaining two pions from the scattering amplitude with a pion, or 
an 7)', propagator. Figure |3| shows the topologies of the diagrams that result. They represent 
the different ways in which flavor can flow in the two step process of pion scattering followed 
by pion propagation. It is only because the coefficient of the chiral log is determined by this 
factorized limit that the flavor flow is well defined. 

Having understood the meaning of the diagrams in Fig. [| I now analyze them in turn. 
The results of this analysis are summarized in Table [I]. The conclusion is that, in the full 
theory, only Fig. |3](a) gives leading order chiral logs in wavefunction renormalization. This 
means that it alone must correspond to the chiral perturbation theory diagram, Fig. Mb). 
In the quenched approximation, however, this diagram is absent. Furthermore, if the quarks 
are degenerate, none of the other diagrams contribute leading order quenched chiral logs. 
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(a) (b) (c) 




(d) (e) (f) 




Figure 3: Quark diagrams corresponding to Fig. 1(b). 



Thus wavefunction renormalization does not contain leading order chiral logs in the quenched 
approximation. 

I begin by analyzing the diagrams in full QCD. It turns out that Figs. |3|(d)-(i) do not 
contribute leading order chiral logs. This is because they all contain one or more mesons 
annihilating into gluons. This can only occur if the meson is the 7/ , because at leading order 
the gluonic amplitude is flavor symmetric. But in full QCD the rj' is massive and its loops do 
not give chiral logs. The only possible contributions are thus from Figs. §(a)-(c). As I now 
explain, however, Figs. |3](b) and (c) give rise only to non-leading chiral logs proportional to 
m^h^m^). 

The extra power of m 2 comes from the four-pion vertex. The two types of diagram that 
contribute to pion scattering are shown in Fig. f|: diagram (a) is required for Fig. |3|(a), 
while diagram (b) is required for Figs. |3](b) and (c). I claim that the amplitude from (a) is 
0(p 2 ), while that from (b) is 0(p A ). (Here I am using p 2 loosely to refer to both p 2 and m 2 .) 
This is true both for the vertex coming from the derivative term in the chiral Lagrangian 
(which is that required for wave function renormalization) and for the vertex from the mass 
term. 

For a theory with four or more flavors this can be seen as follows. One can choose 
the flavors of the external pions such that only Fig. |](b) contributes, e.g. ud, du, sc 
and cs. For these flavors the vertices from the leading order chiral Lagrangian give no 
contribution, because the flavor trace vanishes. Only terms with two flavor traces contribute, 
e.g. Tr^Xc^Xt) 2 , but these contributions are of 0(p A ). For general external flavors, there 
will be contributions to the scattering amplitude from both Figs. (|(a) an d (b). Those from 
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Diagram Full Theory (N > 4) Quenched Approximation 



(a) Leading Absent (no loop) 

(b) Non-leading Non-leading 

(c) Non-leading Absent (no loop) 
(d),(e) Absent (no 7/) Proportional to 5m q 

(f) Absent (no 7/) Non-leading 

(g),(h),(i) Absent (no rf) Absent (no loop) 



Table 1: Contributions of the diagrams of Fig. ^| to chiral logarithms in wavefunction 
renormalization. "Leading" indicates a term proportional to m1ln(m n ), "non-leading" a 
term proportional to vn\ ln(m n ). 
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the disconnected diagram (b) will, however, be of 0(p A ), because the quark diagram is not 
affected if we simply change the flavor labels.^ Thus the only contribution of 0(p 2 ) can be 
from Fig. |](a). That this diagram does contribute at 0(p 2 ) can be seen by choosing flavors 
such that it is the only contribution, e.g. ud, ds, sc and cu. The 0(p 2 ) terms in the chiral 
Lagrangian do give a vertex with these flavors, because they can be connected together in a 
single flavor trace. 

This argument requires only that the leading order interactions are described by the chiral 
Lagrangian. Given the fundamental assumption of this paper, namely that the Lagrangian 
for quenched interactions has the same form as in the full theory, the argument of the previous 
paragraph applies also to the quenched theory. In fact, it applies also if the external particles 
are flavor singlets. Since the background configurations in the quenched approximation are 
by construction independent of the number of flavors of the full theory one starts with, the 
argument also applies to quenched QCD. 

An interesting question is whether Figs. |3](b) and (c) gives rise to leading order chiral 
logs in full QCD. Since QCD has only three light flavors, the argument given above does 
not work: none of the possible external flavors selects Fig. f|(b) alone. This means that the 
question cannot be asked in terms of correlation functions of operators, but only in terms 
of individual contractions. Thus one must have a non-perturbative regulator, such as the 
lattice. Although I do not know of a proof, I think it likely that the N > 4 result carries 
over to N = 3, for the following reasons. First, with staggered fermions, one can show, 
given certain assumptions, that Fig. |](b) is 0(p 4 ) irrespective of the number of dynamical 
quarks PU| . Second, the contribution of wavefunction renormalization in a theory with N 



degenerate quarks, Eq. ([[7]), comes with an explicit factor of N for N > 2. This suggests 
that Fig. |](a) alone is contributing, because there are N flavors of quark running in the 
loop. This is not conclusive, however, because there could be hidden N dependence in / and 
fj,. Finally, the claim that only Fig. g(a) contributes is consistent with the N c — > oo limit, in 
which the both the vertex from the leading order chiral Lagrangian and that from Fig. ^(a) 
are of 0(1/N C ), while those from Figs. |3](b) and (c) are of 0(1/N 2 ). 

I now turn to the analysis of the quark diagrams in the quenched approximation. Of the 
first three diagrams, only Fig. 0(b) is present in the quenched approximation, as the other 
two involve quark loops. As I have just argued, however, Figs. |3|(b) gives non-leading chiral 
logs not only in the full theory (with N > 4), but also in the quenched approximation. Thus 
these three diagrams do not contribute leading order chiral logs to quenched f n . 

The remaining diagrams involve rf propagation. Of these, Figs, ^(g)-(i) require internal 
quark loops, and so are not present in the quenched approximation. Figures |](d)-(f), how- 
ever, do give potential contributions of a type unrelated to those in full QCD. Fortunately, 
if the quarks are degenerate, this contribution vanishes at leading order. This is because, as 
explained in Section |2|, the r( does not appear in the four-pion vertex containing derivatives. 
(Recall that for wavefunction renormalization the vertex must contain derivatives.) For de- 
generate quarks, one can use the usual basis of 7r , r\ and r/, and since the rj' does not couple, 

9 One might worry about the validity of this argument for non-degenerate quarks, for then changing the 
flavor labels also means changing the quark masses. This is not a problem, however, because there is always 
a theory in which, for external flavors chosen so that only the disconnected diagram contributes, the masses 
of the quarks are the same as those in the diagram under consideration. In all such theories the disconnected 
components are 0(p 4 ). 



13 



Figs. |3|(d)-(f) simply vanish. There will, however, be contributions at higher order chiral in 
perturbation theory. 

In fact, Fig. 0(f) does not contribute leading order chiral logs to f n even for non- 
degenerate quarks. This is because only the Vi term in the full chiral Lagrangian, Eq. ([TO]) . 
has the correct flavor structure to correspond to the vertex in Fig. 0(f). Since this term has 
no derivatives, however, it cannot contribute to wavefunction renormalization irrespective of 
the quark masses. 

This is not true of Figs. |3|(d) and (e). For non-degenerate quarks, the neutral eigenstates 
have flavors ss, dd and uu, and each of these individually does couple to the derivative term 
in the Lagrangian. It is simple to see that the contribution to wavefunction renormalization 
vanishes for degenerate quarks because of a cancellation between the two diagrams. For 
non-degenerate quarks, this cancellation does not occur because the neutral eigenstates are 
different. Thus there are extra contributions proportional to the difference in mass of the two 
quarks making up the pion. For example, this extra contribution is proportional to — m u 



for the physical f n . These extra terms are discussed in Ref. and I do not consider them 
further here. 

For mass renormalization Figs. ||](d-f) all contribute even if the quarks are degenerate, 
as discussed further in Section |. 

Finally, I turn to vertex renormalization, Fig. |l](c). The quark diagrams are in one-to- 
one correspondence with those in Fig. |3], and the analysis is almost identical. The result is 
that the summary of Table [3] applies both to vertex and wavefunction renormalization. In 
Fig. I show the diagrams analogous to Figs. 0(a)- (f). I do not show those corresponding 
to Figs. 0(g)- (i), since they contribute neither to full nor to quenched theories. Fig. |(a) 
is the only diagram contributing leading order chiral logs in the full theory, but is absent 
in the quenched approximation. The three pion vertex in Figs. |5|(b) and (c) is absent in 
leading order because it requires two flavor traces, whereas the axial current contains a single 
flavor trace (Eq. (0)). Thus these diagrams contribute only non- leading chiral logs. This 
is true also in the quenched theory, for which only diagram (b) is present. The remaining 
three diagrams require an 7/ in the loop, and are thus absent in the full theory. They are 
present in general in the quenched theory, but vanish for degenerate quark masses because, 
as explained in Section |2], the 7/ does not couple to the leading order representation of the 
axial current. 

In summary, the leading order chiral logs in decay constants, which are present in full 
QCD, are absent in the quenched approximation. There are additional chiral logs from 
7/ loops, which depend on parameters other than / and /1, but these are only present for 



non-degenerate quarks. These extra logs are discussed in Ref. [O 



4 CHIRAL LOGARITHMS IN B- PARAMETERS 

4.1 Bk with dynamical quarks 

To set the stage for the calculation of the chiral logs in Bk in the quenched approximation, 
I review the calculation in theories with dynamical quarks. The discussion will apply to 
theories with a d quark, an s quark, and any number of other light quarks. 
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Figure 5: Quark diagrams corresponding to Fig. 1(c). 



The kaon S-parameter is denned as a ratio of matrix elements 



M K 

M mr . 



(K°\O k \K°) = B K M vac , 

[s a 7 M (l+7 5 )da] [(s67m(!+75)4)] , 

|(^i[^7 M (l+7 5 )4]|0}(0|[(s 6 7 M (l+75)4)]|^°) 



(25) 
(26) 
(27) 



where a and 6 are color indices. The subscript "vac" indicates that this is the result for M.k 
in vacuum saturation approximation. The definition can be written in the more familiar 
form 

M K = fB K f 2 K m K . (28) 



The calculation of chiral logs is clearer using the definition in Eq. (|25|). What is required is 
the logs in both M. K and M. vac . 

I begin with the latter. The matrix element Ai vac is the product of two matrix elements 
of the axial current, the corrections to which are discussed in the previous Section. The 
result can be read off from Eqs. ® and flTBl) 



f 2 Pw ■ p K {\ + 2X(Fig.lb) + 2X(Fig.lc) + 0{m 2 K )) . 



(29) 



The 0(m K ) term comes from non-leading operators, and does not contain chiral logs. As 
in the previous Section, I will drop such terms henceforth. In the formula for X (Eq. |T 



one should use the flavor matrix appropriate to the kaon, = ■^(T e — iT 7 ). The momenta 
are both incoming, so that ■ px = ~Pk = 171 ki where is the physical mass including 
one-loop corrections. The corrections in Eq. ( p9[) are precisely those required to change f 2 
to the one-loop corrected f K 

-Mvac = ffWx ■ (30) 

The chiral logarithmic corrections to M.k in full QCD were first calculated in Ref. 0, and 
related results are given in Ref. [ffijfl . I have checked these calculations, and in the following 
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_ JL_ "0" -Q- SL 

(a) (b) (c) (d) (e) 

Figure 6: Chiral perturbation theory diagrams for M.k- The pair of squares represents the operator 
while the circle is the vertex from the chiral Lagrangian. 

I present a slight generalization of the results. The AS = 2 operator Ok transforms as part 
of a (27, 1) under SU (3) l x SU (3)#. In chiral perturbation theory it is represented by a sum 
of all operators having the same transformation property, with unknown coefficients. There 
are no such operators at 0(1), and only one at 0(p 2 , m q ) |J 

O k - O x K = ±Bf\Y>d^) ds (Y.d^) ds , (31) 

where the subscripts indicate the ds'th element of the matrix. I have chosen the coefficients 
so that, at tree level, Bk = B, as can be seen by combining the tree level result, 

M K = fBf 2 m K , (32) 

with that for A4 v3jC from Eq. (25). Chiral perturbation theory does not determine the value 



of B, which is another parameter like / and /i. 

Chiral logs in M.k come from the five types of diagrams shown in Fig. |6|. I have 
represented the operator by two squares, each corresponding to a left-handed current, because 
it can be written as 

0\ = fB^-iL^L^-iL^) . (33) 

This factorization, which plays an important role in the calculation, is only true at leading 
order in the chiral expansion. The result from these diagrams is a multiplicative correction 

M K = fBf 2 m 2 K (1 + X(Fig.6(a)) + X(Fig.6(b)) + X(Fig.6(c-e))) . (34) 

Here itik is the one loop corrected kaon mass. Strictly speaking, some of the results from 
Figs. H(c)-(e) are proportional to the tree level value olm 2 K (or other tree-level meson masses 
squared). But since these are correction terms, one is justified in replacing m^(tree) with 
mx(lloop). 

The calculation of Figs. |](a) and (b) is identical to that of Figs. [I] (b) and (c) respectively, 
except for an overall factor of two since there are two left handed currents in 0\ 

X(Fig.6(a)) = 2X(Fig.l(b)) ; X(Fig.6(b)) = 2X(Fig.l(c)) . (35) 



This part of the correction to M.k is the same as the complete correction to A^ vac , Eq. (t29|). 
and so cancels in the ratio which defines Bk- In other words, the corrections from Figs. §(a) 
and (b) are precisely those needed to convert the factor of f 2 in Eq. fl34|) to the one loop 
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value f K , and thus do not affect Bk- This is true for any number of flavors in addition to 
the s and d quarks. 

The only new work is that required to calculate Figs. |6|(c)-(e). The result can be written 
in a form which is valid for any number of flavors 



mi 



X(Fig.6(c-e)) = I 2 {m K ) - El 1 + -f)^iK) [(T a )dd ~ (T a ) ss } , (36) 

a m K 



where I2 is the standard integral 



2 

m 



h(m) = — I G(k, mf (37) 



\ 2 



= -L(m)-^—j (1 + OK/V)). (38) 

This integral results from Figs. |6](c) and (d), both of which have a K° in the loop. In Fig. 
|6|(e) there are only flavor-diagonal pions in the loop. This diagram gives rise to the second 
term in X, in which the sum is over flavor-diagonal pions in the mass eigenstate basis.0 

An important consistency check on this result is that divergent parts can be absorbed 
into B. This requires that there be no quartic divergences, for these would imply that 
M.K oc A 4 , while chiral symmetry dictates that the leading term is proportional to m 2 K . 
All three diagrams do contains quartic divergences, but they cancel in the sum. There are 
quadratic divergences in Eq. (|36| ) because Ii oc A 2 at leading order. These can be absorbed 
into B as long as they do not depend on the quark masses. The factor of m 2 a /m 2 K multiplying 
Ii{m a ) suggests that this may not be true, but in fact the mass dependence cancels because 
of the result 

J2mH(T a ) dd -(T a ) ss } 2 = m 2 K . (39) 

a 

Here again the sum is in the mass eigenstate basis, and only flavor-diagonal states contribute. 
This result holds for any number of flavors. 

Combining Eqs. fl30|), (H), (|35]) and ([§), the corrected B K is 

B K = B { 1 + I 2 {m K ) - £(1 + ^)h(m a ) \{T a ) M - (T a ) s A . (40) 

In QCD (with m u = m^) this gives 

(TYl 2 Tfl 2 \ 

1 + h(m K ) - |(1 + -f )IiK) - |(1 + ^ L )/i(m 7r ) . (41) 
m K m K ) 

In the limit m u = m d — > the chiral logarithmic part is 

B K ^B(l-AL(m K )) , (42) 



'in fact, the sum can be extended to all pions, because the last factor is non-zero only for flavor-diagonal 



pions. The same is true in Eq. (40) below 
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where the quadratic divergences have been absorbed into B. It is also useful to expand about 
m s = md, in terms of 5 = (m s — m^) 2 / \m s +md) 2 

B K = B(l - L(m K )(3 + §5) + 0(5 2 )) , (43) 

a result previously quoted in Ref. Hl5| . Finally, for N degenerate flavors one finds 

B K = B(l + I 2 {m) - 2h(m)) = B(l - 3L(m) + ...). (44) 

Unlike the correction to /„■, this result has no explicit dependence on N. 



4.2 Bk in the quenched approximation 

I now turn to the calculation of the chiral logs in quenched Bk- The conclusion of the 
following analysis is simple: the result for Bk is unchanged from the full theory except that 
(1) there are contributions from rf loops when m s ^ rrid, and (2) in the formula Eq. fl40|) the 
mass eigenstates are ss and dd (with masses m ss = 2fim s and m^d = ^firrid, respectively). 
Thus the quenched result is 

^(quenched) = B (l + h (m K ) - ^^hM - 

+(m a -m d )X(rf)) , (45) 

where X(rj') comes from rf loops. In the limit = m s the result is identical to that in the 
full theory, Eq. (@§. 

To begin the quenched calculation, the operator Ok must be written in terms of pion 
fields. As before, I assume that the form of the operator is the same as for the full theory, Eq. 
(|31"D, except that the constant B can differ. In other words, I assume that the amplitudes for 
the operator to create two and four pions are related in the same way as in the full theory. In 
the continuum, these relations are derived from the chiral Ward identities. With staggered 



fermions, analogous Ward identities hold in both full and quenched theories [|19i |15] . With 
Wilson fermions, the Ward identities are regained in the continuum limit of both full and 
quenched theories Dynamical fermions play no role in the derivations. 

Given this assumption, the calculation proceeds as in the continuum, except for the usual 
two complications. First, one must keep only those diagrams in Fig. |6] which correspond to 
quark diagrams without internal loops. Second, one must determine whether there are extra 
contributions from r/ mesons propagating in the loops. 

For Figs. |5](a) and (b) the analysis from the previous Section can be used, because the 
operator 0\ is a product of two factors of L^. The calculation of Figs |6](a) and (b) is 
identical to that of Figs. |l](b) and (c). This means that, just as in the full theory, any 
chiral logs from Figs. |6](a) and (b) cancel in ratio which defines Bk- In fact, as shown in 
the previous Section, the chiral logs in quenched Jk differ from those in the full theory. In 
particular there are no logs in Jk for m s = rrid, and there are extra logs from r/ loops if 
m s 7^ m d- But these differences are irrelevant for B K - 

Just as in the full theory, then, the logs for Bk come entirely from Figs. |6|(c)-(e). The 
quark diagrams corresponding to these are shown in Fig. [7[ The operator Ok is shown 
as two squares, each representing a left handed current. Each diagram actually represents 
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two contractions, one in which the color indices are contracted within each current, the 
other having the indices contracted between currents. The pairings of quark and antiquark 
propagators into bilinears can be interchanged by a Fierz transformation. I have chosen the 
pairing so that each bilinear couples only to one of the mesons in the loop. This makes it 
easier to see the correspondences with pion diagrams. 

None of the diagrams in Fig. contain internal quark loops. Thus all are present in both 
quenched and full theories. The correspondences between quark and pion diagrams is as 
follows. Figures 0(a) and (b) correspond to Fig. |6|(c), since both have a four pion scattering 
amplitude. Figure 0(c) corresponds to Fig. H(d), since the particle in the loop is a kaon. 
And, finally, Figs. 0(d) and (e) correspond to Fig. |6|(e), in which the mesons in the loop 
have flavors ss and dd. 

Since Figs. 0(a), (b) and (c) are the same in both full and quenched theories, and since 
the particles in the loops are flavor off-diagonal, the result from the chiral perturbation 
theory diagrams Figs. ^|(c) and (d) is identical in full and quenched theories. This is the 
first example in which there are chiral logs in a quenched quantity coming from the same 
diagrams as those in the full theory. In fact, Fig. 0(b) does not contribute at leading order 
because the four pion vertex is 0{p A ). 

The quenched approximation can, however, affect the evaluation of Fig. |](e), because the 
loop contains flavor- diagonal mass eigenstates, which, in general, are different in quenched 
and full QCD. There are also potential r^'-loop contributions in the quenched theory not 
present in the full theory, for example that of Fig. 0(e). The effect of these differences is 
that, in the expression for X(Fig.|6](c — e)), Eq. (|36D , the sum in the quenched theory runs 
over both the pions and the rf , while in full QCD the rj' is excluded. The difference between 
these two sums vanishes, however, for degenerate quarks. This is for two reasons. First, the 
non-singlet pion mass eigenstates are the same in quenched and full theories, so that the 
contribution to the sums from these pions agree. Second, the rf is a mass eigenstate, and 
its extra contribution to quenched sum vanishes, because T ss — T^d = for a flavor singlet. 
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Another way of seeing this is to recall that the rf does not couple to L^, so the extra diagram 
0(e) vanishes. 

It follows from this discussion that Eq. (|45|) describes the leading order chiral logs in 
quenched B K . 



The matrix element Mk picks out the positive parity part of Or, which consists of a sum 
of eight different terms: four "vector- vector" (7^ • 7^) and four "axial-axial" (7^75 • 7^75). 
There are two different color contractions for each of these eight terms, with one or two 
traces over color indices. When doing a numerical calculation one calculates all sixteen terms 
separately, so it is interesting to predict their chiral behavior and finite volume dependence. 
In this subsection I calculate the leading order chiral logarithms for the four Lorentz invariant 
combinations in which the index \x is summed. 

The precise definition of the breakup of Mk into sixteen pieces is this: 

1. Write down all contractions (i.e. products of quark propagators). 

2. Fierz transform the operator such that both external kaons are connected to one or 
other bilinear. All terms will thus have two traces over spinor indices, and either one 
or two traces over color indices. 

3. Break up each of the resulting terms into the eight vector- vector and axial-axial pieces. 

Clearly, this definition applies at the level of contractions. In order to use chiral perturbation 
theory, however, one must have a definition in terms of matrix elements of operators. For 
three or less flavors it is not possible to write down operators whose matrix element give 
each of the sixteen terms separately. As I show shortly, to do this requires at least four 
valence flavors. In quenched QCD, this presents no problem, for one has already decoupled 
the number of valence and dynamical flavors. It does mean that the results derived below 
are not directly applicable to full QCD. They should, however, be as valid for full QCD as 
they are for quenched QCD. 

Consider, then, a theory with N > 4 light flavors: an s and an s', both of mass m s ; a d 
and a d', both of mass m^; and N — 4 other quarks with arbitrary (though small) masses. 
Let the K be the ds pion, and the K' be the corresponding state made out of primed quarks. 
Define Mk and Bk in this theory exactly as in QCD 



4.3 By and B A 




(46) 
(47) 



16 o f2 2 



Then, by comparing contractions, we can write Mk as follows 



M K 
M v 

M A 
M V i 



2(M v + Ma) 
M V i + M V2 
Mai + M A 2 



(48) 
(49) 
(50) 
(51) 
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M V 2 = m[s' alfM d' a ] \s bl MK) (52) 
M Al = (IC\ [s' aW d' b ] [s bW d a ]\K) (53) 
M A2 = (K 7 \[s' alfll5 d' a ) [s bl , l5 d b )\K) . (54) 

The point of the extra flavors is to restrict contractions; each of the matrix elements in the 
last four lines has only one contraction. The factor of two on the first line appears because 
the two bilinears in Aix are identical, and so there are twice as many contractions as in M.y 
and M.A- The subscripts 1 and 2 indicate the number of traces over color indices. 
The corresponding .B-parameters are defined by 

Mvi = HWkBvx , M V 2 = lf 2 K m 2 K Bv2 , (55) 

and identical equations with V — > A. These definitions are such that 

B K = B V + B A = B V i + B V2 + B A1 + B A2 . (56) 

It is straightforward to show that in vacuum saturation approximation B A \ = 0.25, B A2 = 
0.75, and thus B A = 1, while By = Byi = By 2 = 0. The numerical results, shown in Section 
0, do not, in fact, agree with these predictions. 

I calculate the chiral logs in these i?-parameters in the full theory, and then argue that 
the results apply with only minor changes in the quenched approximation. The first step 
is to transcribe the problem into the language of chiral perturbation theory. As explained 
in section |^, the chiral Lagrangian itself takes the same form for any number of flavors. To 
transcribe the operators, however, is more complicated. It turns out that, at the order I am 
working, one needs the following set 

1 = f(W^) d , s ,(Zd^) ds (57) 

2 = / 4 (S^St) d , s (S^St) ds , (58) 

3 = -2/i 2 / 4 £^4 s , (59) 

04 = f\d^ dls ){d^ ds ,) (60) 

O, = / 4 (£^£tW(£t^£) ds . (61) 

The overall factors are chosen for later convenience. 

To transcribe the quark operators, they must be written in terms of representations of 
the SU(N)i x SU(N)r chiral group. The first step is to rewrite the operators in terms of 
left- and right-handed bilinears 

[V'V\i = l(lL'L]i + [L'R]i) + O- , [A'A]i = i([L'L], - [L'%) + 0~ , (62) 

where i — 1, 2 labels the number of color traces. The negative parity operators 0~ and 0'~ 

do not contribute to the matrix elements under study and can be ignored. I am using the 
compact notation 

[W]i = KlA\ %lM > (63) 

[A'A]! = [s / 7m7b4] \sbl»l$da] , (64) 

[L'L]i = [<7m(1+75)<] [s 6 7„(1+75)4] , (65) 

[L'R]i = Kl^+lM [s&7m(1-75K] , (66) 
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with similar definitions for operators with two color traces. 

The left-handed bilinears L and 11 transform as ((N 2 — 1),1), i.e. as adjoints under 
SU(N) L and singlets under SU(N) R . Under SU(N) L , [L'L] lies in the symmetric product 
of two adjoints, 

((N 2 -l) ® (N 2 -l)) symm = 1 + (N 2 -l) + S + A . (67) 

The representations S and A are both traceless, and are respectively symmetric and anti- 
symmetric under interchange of quark indices. They have dimensions (N — 1)N 2 (N + 3)/4 
and (N — 3)N 2 (N + l)/4 respectively. Since 11 and L have no flavors in common, their 
product must be a combination of S and A. The combinations belonging to one or other 
representation are 

\L'L\ 2 ±\L'L\ X = [s' al ^ 5 )d' a ] [s bl ,(l^ 5 )d b ]±[s' al ,(l^ 5 )d a ] [stf^l+yM e I ^ . (68) 

In chiral perturbation theory the lowest dimension operators transforming in these represen- 
tations are 

(<S, 1) : O s = 1 + 2 ; 
{A, 1) : O a = 1 - 2 . 

Thus the transcription of the continuum operators is 

[L'L^ + iL'Ll = a s O s + 0(p 4 ) , (69) 
[L'L] 2 — [L'L]i = a A A + 0(p 4 ), (70) 

where as and a A are new, unknown coefficients. One can estimate their values by taking 
matrix elements between a K and a K', using the vacuum saturation approximation for the 
l.h.s., and tree level chiral perturbation theory for the r.h.s. The resulting estimates are 
ctg — g and ol a — fr. The same results are obtained using external d s and s'd kaons. 

The [L'R] operators transform as ((N 2 — 1), (N 2 — 1)). They can thus be written as linear 
combinations of the three operators O3-5, which also transform in this way 

[L'Rji = fc0 3 + lt O A + 5 % O b + 0(p 4 ) . (71) 

There are six unknown coefficients, three for each color contraction. To estimate their size 
I first take the matrix element between a K and a K' . At tree level, only O5 has a non- 
vanishing matrix element. Using vacuum saturation for the matrix element of [L'R]i yields 
the estimates 5i — |, 62 — 1. 

To estimate the $ I take the matrix element between d's and s'd kaons. In chiral 
perturbation theory one finds, at tree level, 

(K(d's)\[L'R} t \K(s'd)) = -mfi 2 f 2 + 0( lt m 2 K ) . (72) 

To use the vacuum saturation approximation I first Fierz transform the operators, e.g. 

[L'R], = -2^(1-75)4] fc,(l+7 5 K] , (73) 
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Operator 
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c 3 


c 4 


c 5 


2[V'V] 1 




a + 
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7i 


5i 


2[V'V) 2 


a + 


a_ 




72 


5 2 


2[A!A\ 1 


a_ 


a + 


-A 


-7i 


-5 


2 [AM] 2 


a + 


a_ 


-A 


-72 


-5; 



Table 2: Decomposition of quark operators into chiral operators J2j c jOj- 
and then use 

(0\s b ^d' b \K(^d)) = ifiV2f K , (74) 

to give 

(K(d's)\[L'R] 1 \K(s'd)) vac = -4/i 2 / 2 , (75) 
(K(d!s)\[L'R} 2 \K(s'd)) yac = -|/i 2 / 2 • (76) 

The resulting estimates are Pi = 1 and /3 2 = §■ This method gives no estimate for the 7j. 

In summary, each quark operator can be written as J2j c j@j with corrections of 0(p 4 ). 
The coefficients are collected in Table where I use the definitions 

«± = 5 («s ± cka) • (77) 

In vacuum saturation approximation these are a + — 1 and cn_ = |. 

All four operators in Eqs. (|5l"|)-(|5"4]) contain "left-left" and "left-right" parts. Equation 
( |72|) shows that the matrix elements of "left-right" operators are generically of 0(1), unlike 
those of "left-left" operators which are of 0(p 2 ). Nevertheless, at tree level, and for the 
external flavors K and K', all four operators have matrix elements of 0(p 2 ). The choice 
of external flavors has removed the 0(1) contributions. As we will see, however, these 
contributions do have an effect when loops are included. 

The calculation of the chiral logs parallels that for Bk- The same diagrams contribute 
and the general result is 

(IC^lcjOjlK ) = ( Cl -c 5 )2m 2 K f K 

+ c 2 f \2m 2 K h{m K ) - (m 2 K + m 2 s )Ji(m ss ) - (m 2 K + m^)/i(m dd )) 
+ c 3 2/i 2 / 2 (h(m ss ) + h(m dd ) - 2h(m K ) - 2I 2 (m K )) 
+ c 4 f (Am 2 K h(m K ) - 2m 2 K I 2 {m K ) 

-{m 2 K + m 2 ss )h(m ss ) - (m 2 K + m 2 M )h{m dd )) . (78) 

I have simplified the result using m s r = m s and m d > = m d . Only the operators 0\ and 
5 contribute to Figs. ||(a) and (b), and the corrections simply convert f 2 into the 1-loop 
result f\. The general expression for f K at 1-loop is given in Section [| These are the only 
contributions of 0\ and O^. The remaining three operators contribute only to Figs. |](c)-(e). 
The pions in the loops in these diagrams always consist of a primed quark and an unprimed 
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antiquark (or vica- versa). These states are automatically mass eigenstates, for any number 
of flavors, so the result is valid for all N > A. 

There are no quartic divergences in Eq. (|78|). The only quadratic divergences at this order 
are proportional to c 2 vr? K A 2 . These can be absorbed into the coefficient c\, and correspond 
to mixing between 0\ and 2 - The quadratic divergences proportional to C3 cancel; this 
is essential, because they are not proportional to m 2 K , and could not be absorbed by a 
redefinition of the coefficients. The quadratic divergences proportional to C4 cancel because 
2771^ = m 2 ss + m 2 dd . This is not required by chiral symmetry, however, as a divergence could 
have been absorbed into C5. Finally, there are also quadratic divergences hidden in f^, which 
can be absorbed into /, as discussed previously. 

1 now collect the results for the S-parameters, using the general formula Eq. (|7H|), the 
decompositions given in Table [| and the definitions Eqs. (^8|) -(^). I present the results 
for m s = rrid, since it turns out that these also apply in the quenched approximation. The 
generalizations to non-degenerate quarks are easy to obtain from the above. For the vector 
5-parameters I find 

Bvi = |(a_-<y 1 ) + |a + (/ 2 -2/ 1 )-f/3 1 -^-/ 2 -| Tl / sl , (79) 

m K 

B V2 = |( a+ - ( J 2 ) + | a _(/ 2 -2/ 1 )-f/? 2 A-/ 2 -|7 2 / 2 . (80) 

Since mx = rn ss = rridd, the arguments of the Ii are all the same, and I have dropped them. 
The results for Ba are the same except that terms proportional to /3j, ji and Si have the 
opposite sign. A check on these results is that they combine to give the same result for 
B K = B V x + B V2 + Bay + Ba2 as found previously (Eq. (p4|)) 

B K =\a s (\ + h-2I x ) , (81) 

providing one identifies B = 

The corrections to By and Ba include a term of the form l2/m 2 K which is proportional 
to ln(m#). This is enhanced by a factor of fi 2 /m 2 K over the usual chiral logs. The general 
possibility of enhanced chiral logs was noticed long ago by Pagels and Langacker M. The 
enhanced log here is the remnant of the fact that O3 is of 0(1) in the chiral expansion, 
rather than of 0(p 2 ). For small enough mx the enhanced log dominates over the "leading 
order" terms, and it diverges in the chiral limit. The vector- vector and axial- axial matrix 
elements do not themselves diverge, but contain a chiral log proportional to m\\n{rriK), 
rather than the usual m^-ln(m^-). The unknown scale in the log can be absorbed into Si, 
and corresponds to mixing between O3 and O5. 

These results can be adapted to the quenched approximation with the by now standard 
procedure. The reasoning is identical to that for B K , because the quark diagrams are the 
same. The only diagrams containing internal quark loops are those that change f 2 to f\, 
and these cancel in B— parameters. The other diagrams, those of Fig. 0, are included in the 
quenched approximation. There are only two possible differences between quenched and full 
theories. The first comes from the fact that the mass eigenstates are not the same for flavor 
diagonal states. This is not relevant here, however, since all pions in the loops are flavor 
off-diagonal. The second comes from rj' loops, which do not contribute to the full theory, 
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but might be present in the quenched approximation. The quark diagram is that shown 
in Fig. 0(e). As for Bk, however, this diagram does not contribute if m s = rrid- For 0\, 
2 and 5 , this is because the rf does not couple to these operators. For O s and O4, it is 
because the rj', although it couples to the operators, cannot appear in a loop like that in 
Fig. 0(a) because the flavor indices are wrong. The rf does appear in two loop diagrams. 
The conclusion is that, for m s = ma, the results of Eqs. ([78|)-(|5D|) apply without change to 
the quenched approximation. For m s 7^ m^, however, there will be additional contributions 
from Fig. 0(e). 

There is one complication in the numerical calculation of the B— parameters which arises 
if one uses staggered fermions. I describe this briefly here; for more details see, for example, 



Refs. [pL9|, |15| . Each lattice staggered fermion represents four continuum flavors, so the quark 
fields have an additional flavor index. The bilinears that are actually used have these extra 
indices contracted as in 

sin d — > Si7,A + s 2 7^2 - s 3 ^d 3 - 347^4 , (82) 

where the color indices are implicit. This particular flavor is that of the lattice Goldstone 
pion, and in numerical calculations the external kaons have this flavor. The analysis given 
above is valid for any number of flavors, and thus applies directly to these generalized 
operators. One simply has to break each matrix element up into parts, each of which is 
like those in Eqs. (^TJ)-(^). Keeping track of flavor factors, one finds that the form of the 
final result is unchanged. The only difference is that most of the pions in the loops have 
flavors such as s^d^, and thus are not Goldstone pions. In fact, it is straightforward to show 
that only one out of sixteen pions is a Goldstone pion. This is important at finite lattice 
spacing, because non-Goldstone pions are slightly more massive than the Goldstone pion, 
and one must use the appropriate mass in the loop integrals. In addition, the vertices of the 
non-Goldstone pions will differ from those of the Goldstone pions by terms of 0(a). The 
effect of this is to make the replacement 

h{m K )/m 2 K — ► ±I 2 {m G )/m G + ^R N Gh{m NG )/m 2 NG , (83) 

where m G and m^ G are > respectively, the mass of the Goldstone and non-Goldstone pions, 
and Rng is the ratio of non-Goldstone to Goldstone coefficients. (I am assuming here that 
all the non-Goldstone pions are degenerate, which, numerically, is a good approximation.) 
Similar replacements should be made in all the loop integrals in Eqs. (]79f) and flSOl). 



5 7/ LOOPS 

The major new feature of quenched chiral perturbation theory is the appearance of the light 
rf in loops. Thus far I have considered only quantities which do not have rf contributions at 
one-loop order. Clearly it would be preferable to understand such contributions, since this 
would increase the number of quantities whose finite volume dependence could be predicted. 
In addition, it is potentially dangerous to leave an unsettled problem in one part of the 
theory, as it may feed back at higher loop order into the quantities I have studied. 

I consider the effects of rf loops in theory with iV degenerate valence quarks. Using 
degenerate quarks simplifies the discussion in two ways. First, the rf can be chosen as one of 
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(a) (b) (c) 

Fi gure 8: Diagrams contributing to mass renormalization in the quenched approximation. The circle 
represents the vertex from the mass term in the chiral Lagrangian, the cross an insertion of the gluonic 
amplitude C n '- Lines with a cross are rf propagators; the others can be either pion or rf propagators. 

the mass eigenstates, which allows the effects of rf loops to be clearly distinguished. Second, 
as discussed above, the rf only couples to the pions through the mass term in the chiral 
Lagrangian, C m . This means, in particular, that rf loops only affect mass renormalization, 
and I focus on this in the following discussion. 

In the quenched approximation, the only diagrams which renormalize the pion mass at 
leading order are those of Figs. |3](d)-(f). Note that only the rf appears in these loops, since 
non-singlet states cannot annihilate into gluons. As discussed at the end of Section |^, the 
contribution from Fig. [|(f), which arises from the vertex V^O)?/ Tr(M?/7r 2 ), can be rotated 
away by a field redefinition. Its effect appears in a change in the value of the constant A in 
C n ' (Eq. Q14D). Thus I need only consider Figs. ||(d) and (e). These are represented in chiral 
perturbation theory by Fig. §(a), where the cross is the gluon annihilation vertex, C^. The 
renormalized mass at one-loop from Fig ||(a) is 

ml = 2/xm(l +X(Fig.8a)) (84) 
X(Fig.8a) = [G(k,m rf ) 2 (Am 2 Q + (A- l)k 2 ) (85) 

-(A- 1)) J 2 K0 + 1(4 - l)IiK0 (86) 




The mass in the loop is m' which equals m w at tree level. This result agrees with Ref. 
The k 2 term in the integrand leads to a correction of the usual type: it is proportional to 
m 2 t ln(m^'), and vanishes in the chiral limit. The m\ term, by contrast, leads to a correction 
proportional to m^Xivim^ which diverges in the chiral limit. This enhanced chiral log is 
similar those appearing in By and The additional dimensionful parameter m violates 
the standard power counting rules of chiral perturbation theory, according to which each 
loop comes with a factor of m 2 , possibly multiplied by logarithms. For small enough m n , 
the one-loop correction proportional to mgk^m^/) is not small, and one must sum all such 
terms. 

In practice, as stressed in Ref. [iT|, this summation may not be numerical important. 



This is because in present lattice calculations the pions are not very light, m n > 400 MeV. 
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To estimate the size of the correction, I extract the leading logarithm 




X(Fig.8a) = -6\n[ ) . (87) 
Here I have introduced the parameter 



Ami 

~ 0.2 



N(4nfy 

where the estimate uses N = 3, A = 1, and mo = 900 MeV. Taking A »s m p , the magnitude 
of the logarithm is roughly unity for m n = 400 MeV. Thus the correction increases m\ by 
only 20%, and resummation is not required. Since the logarithm is not large, however, it 
does not dominate over constants and power corrections, and the full expression of Eq. (j86|) 
should be used. 

Bernard and Golterman have begun a program of calculations, in which they include 
rf loops only at one- loop order, for the general case of non-degenerate quarks [Tl||. Their 
predictions can be compared to numerical results for various quantities. My interest in this 
section is not phenomenological, but rather to understand the apparent divergence of the 
correction in the chiral limit. To address this question one must resum the leading logarithms 
of the form m 2 / 1 \n(m rj i) n / f 2n . This is relatively straightforward, because the leading logs 
come only from diagrams in which £ v t appears on a propagator both ends of which come 
from the same vertex. An example is shown in Fig. ||(b). Diagrams involving rf propagators 
connecting two or more vertices have fewer logs. For example, the diagram of Fig. H(c) gives 
a correction proportional to m^m 2 ,// 6 , with no logs at all. 

The bubbles exemplified by Figs. |^(a) and (b) sum up to an exponential 

ml = 2/jmexp (X(Fig.8a)) , (89) 
the leading logarithmic part of which is 

/ A 2 \ S 

ml = 2/im —5- . (90) 



It is straightforward to show that the same exponential factor corrects any vertex coming 
from C m , so that the entire mass term is replaced by 




C m — 2 ) t-"m j (91) 



with the proviso that no 77' fields coming from this vertex be joined together by an insertion 
of £ v i. One corollary of this result is that the rf mass itself is corrected in the same way 
as TOtt, so that they remain equal. This means that Eq. ( |9"0| ) should be interpreted as a 
self-consistent equation for the common pion and rf mass, with m v > on the r.h.s. replaced 
by m n . Solving this equation yields the final result 



s 

ml 



(A 2 ) 1 "* . (92) 
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(a) 




Figure 9: An example of a diagram corresponding to self-consistent leading-log evaluation of m T in the 
quenched approximation: (a) in chiral perturbation theory; (b) in terms of quarks. Notation as in Fig. 8. 

Solving Eq. (R^) self-consistently corresponds to summing all the diagrams of the form Fig. 
|9|(a), in which each internal rj' propagator is decorated with any number of bubble sums. In 
terms of quarks, these "cactus" diagrams look like Fig. |9](b). 

The result Eq. (|92]) seems very odd at first sight. For one thing, m 2 has a non- analytic 
dependence on the quark mass. For another, it appears to contradict the strong numerical 
evidence that m 2 oc m in quenched as well as full theories To investigate further, I 

substitute Eq. (|92]) back into Eq. (fH|) to give (using also Eq. (|13|)) 

f 2 i -J- 

C = —T ( 2 V m )~ ( a2 ) 1+5 (cos(0 o ) Tr + + isin(0 o ) Tr [£-£ f ]) . (93) 

This says that, at leading- log order, the effect of rf loops amounts to a change in the coeffi- 
cient multiplying the mass term. The quenched chiral Lagrangian looks the same as that in 
the full theory Q as long as one uses a new mass m! 

( A 2 \^ 

m = m \ . (94) 

\2fimJ K ' 

Looked at this way, the result is less strange. All the relationships between amplitudes 
involving different numbers of pions remain valid, e.g. the Weinberg scattering length for- 
mulae. These relations involve only the physical parameters f n and m n , and not the quark 
masses. In particular, there is no contradiction with exact lattice Ward identities that can 



be derived with staggered fermions [19, 20 



Associated with the enhanced chiral logarithms is an "enhanced" finite volume depen- 
dence. The precise form of this dependence is discussed in section 7. The important point 
for present purposes is that it too is as a function of physical quantities alone, and not of 
m (or m'). Thus one cannot use the finite volume behavior to observe the non-analytic 
dependence on m. 

Is there any way of observing this dependence? The remaining hope rests with the 
numerical lattice results for the quark mass dependence of m 2 and the condensate (ipip)- 

11 Of course, in the application of this Lagrangian, one must follow the rules discussed in the sections 3 
and 4 for discarding diagrams containing internal quark loops. 
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The predicted dependence of m\ on m is given in Eq. fl92"l) . That of the condensate may be 
deduced from Eq. (or from the formula of Gell-Mann, Oakes and Renner) 



{ipip) (xm^ 6 . (95) 
These predictions are in apparent conflict with numerical simulations which find that m\ oc 



m q and (ipi/)) = (V'V ; )o + 0(m q ), where m q is the bare lattice quark mass p4|. The results 
are probably good enough to rule out 5 ~ 0.2, assuming that we identify m q with m. 

In full QCD this identification is needed for the following reason. The bare quark mass 
m q is part of a scalar source s. Functional derivatives with respect to s yield correlators 
involving scalar bilinears in the presence of the source. Setting the source to zero we obtain 
the correlators in the massless theory, the condensate being the simplest example. We expect 
these correlators to exist in a regulated theory. For this to be so, the chiral Lagrangian must 
contain only positive integral powers of s. This means that, up to an overall constant, and 
up to corrections of 0(m 2 ), we should identify m q with m. 

In quenched QCD, by contrast, one does not have a Lagrangian framework, and I see no 
reason why the connection between m and m q should be analytic. I propose that m q should 
be identified with ml (up to an overall factor), rather than with m. This would imply the 
relationship 

m q oc m 1+s . (96) 

As far as I can tell, once this identification is made there are no observable effects of the 
non-analytic dependence on m. 



6 FINITE VOLUME DEPENDENCE 

The discussion thus far concerns quantities calculated in infinite volume. Numerical simu- 
lations, by contrast, are done in finite boxes. These have finite extent in both space and in 
Euclidean "time". I denote the length of the lattice in the spatial directions by Lj, i = 1—3, 
and the length in the time direction by L 4 . The latter is related to the temperature by 
T = I/L4. In the following discussion, I assume that L 4 >> Lj, so that the dominant finite 
size effect is due to the spatial volume. 

Gasser and Leutwyler have shown that the same loop graphs which give rise to the leading 
chiral logarithms also allow one to predict the leading finite volume dependence ||. They 
have presented results for f n and in a theory with N degenerate flavors. It is simple to 
extend their method to any quantity for which the chiral logs are known. The only subtlety 
concerns the applicability of their method to the quenched approximation. 

Their method depends upon the result that the volume dependence in chiral perturbation 
theory enters only through the pion propagators. The chiral Lagrangian itself does not 
depend on the volume, and so the vertices in chiral perturbation theory are unchanged. This 
is also true of external sources, such as that representing Ok- Thus the constants /, /1, and 
B are volume independent. The volume dependence of the pion propagator enters through 
the boundary conditions, which are periodic in all four Euclidean directions. This means 
that the integral over momenta is replaced by a sum, in which is an integer multiple of 
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27r/L Al . In position space, the propagator is a sum over periodic images 

G L (x, m) = G(x, m) + ^2 G(x + n^L^, m) , 



(97) 



where = (1, 0, 0, 0), (—1, 0, 0, 0), (0, 1, 0, 0), etc. Volume dependence comes from replacing 
G with Gl in loop integrals such as l\ and li- 

The volume independence of the chiral Lagrangian is clearly crucial to the method. 



that chiral Ward identities remain valid in finite volume, (b) causality, (c) that the pion 
propagator at fixed distance should approach the infinite volume form as the volume increases 
to infinity, and (d) the permutation symmetry between the four directions /i = 1 — 4. The 
first three ingredients establish the existence of a finite volume chiral Lagrangian, and that 
the coefficients of its leading terms are the same as at infinite volume. Volume dependence 
suppressed by powers of 1/L 2 is allowed, so there can be terms such as J2i Tr^E^E^/L^. 
Such terms are ruled out by the fourth ingredient, permutation symmetry. This symmetry 
implies that, if there is any dependence on the spatial lengths Lj, there must also be a 
dependence on L 4 , and thus on T. But this contradicts a standard field theoretic result, 
namely that the temperature enters the functional integral representation of the partition 
function only through the boundary conditions. Thus permutation symmetry rules out Li 
dependence in the vertices at any order in the chiral expansion. For the symmetry to be 
present, the boundary conditions must be the same in all four directions. In particular, 
since the finite temperature functional integral has antiperiodic boundary conditions on 
fermions in the time direction, the boundary conditions must also be antiperiodic in the 
spatial directions. 

In the quenched approximation only the first and third ingredients are present. Concern- 
ing the former, one can derive the Ward identities on the lattice independent of the quenched 
approximation [18|, [L9| . As for the latter, it is reasonable to assume that the boundaries be- 
come unimportant as they recede to infinity. But the lack of dynamical quarks means that 
one cannot even formulate the second ingredient, causality. Causality requires that the 
anticommutator of fermion fields vanish outside the forward lightcone. In the quenched ap- 
proximation, one cannot construct the transfer matrix and obtain an operator formulation in 
which quarks are represented by anticommuting operators. A less formal way of seeing this 
is to note that, without quark loops, there is no Pauli exclusion principle, so that, in some 
sense, fermion fields no longer anticommute. The lack of quark loops also means that the 
fourth ingredient, permutation symmetry, is not relevant. The requirement of antiperiodic 
boundary conditions applies to internal quark loops, for it is these loops which build up 
the contribution of states containing quarks to the partition function. Since these loops are 
absent, one is free to choose the boundary conditions on the valence quarks as one pleases. 

The absence of two of the four ingredients means that that one cannot rigorously argue 
for the volume independence of the chiral Lagrangian in the quenched approximation. In 
fact, as already discussed, even the existence of such a Lagrangian is an assumption. I will 
simply extend this assumption by taking the quenched chiral Lagrangian to have volume 
independent coefficients. I think this is reasonable, and in the following I give a qualitative 
argument in its support. Ultimately, however, it must be tested by numerical simulations. 

12 In addition, one must work at volumes large enough that chiral symmetry is not restored. 



The arguments for this are presented 
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The physical origin of the finite size effects from chiral loops is that particles are sur- 
rounded by a cloud of virtual pions. The pions in this cloud can propagate "around the 
world" due to the periodicity, the propagator falling roughly like exp(— m n Li). If there were 
no interactions, then there would be no pion cloud, and the only effect of finite volume 
would be to restrict the momenta to discrete values. In particular, masses would remain 
unchanged. Luscher has given a general analysis of finite volume effects arising from polar- 
ization clouds on the masses of stable particles ]?]]. For the pion, his analysis implies that 
the leading order effect comes from a "tadpole" diagram, similar to that of Fig. |l|(b), except 
that the vertex and propagator are fully dressed. The crucial result is that volume depen- 
dence enters only through the pion propagator, which goes "around the world" in each of the 
six spatial directions. Liischer's general formula simplifies when applied to pions, because 
one can do a chiral expansion of the propagator and vertices.^ The result coincides with 
that obtained from chiral loops, up to exponentially small corrections which are beyond the 
accuracy of Liischer's formula ||. This agreement would not hold if there were any volume 
dependence in the coefficients of the chiral Lagrangian, such as terms of 0(1/L 2 ). This 
provides independent confirmation of the arguments of Ref. |25] . 



Liischer's result can be paraphrased as follows. The pion only "knows" about the finite 
size of the box through loop diagrams, and in particular through the pion propagating in the 
loop. The vertex which emits and absorbs the pion does not itself "know" about the finite 
volume. Such knowledge would require another pion loop, which would bring an additional 
exponential suppression. This physical picture applies also to properties of the pion other 
than its mass, for example /„- and Bk- It also seems reasonable to apply it to the quenched 
approximation. 

The previous discussion has ignored the possibility of finite volume effects due to particles 
being "squeezed" by the box. To avoid such squeezing the box must be larger than diameter 
of the particle. It is not very clear what box size this requires. Numerical calculations 
of various definitions of "wavef unctions" show that they fall approximately exponentially 
in the range 0.5-1.5 fm, with a scale of roughly m-yvF = 750 MeV |26|, |27| . This may 
imply a correction factor from squeezing that falls as exp(— m^pL). As long as < 
wwf, the contribution from squeezing should be suppressed. At large enough distances the 
wavefunction should fall faster than an exponential, so the effect of squeezing should be 
further suppressed. 

Another source of finite volume dependence is loops of heavier mesons, such as the p. 
Such contributions would fall roughly as exp(— m p L). As long as << m p , they should 
be strongly suppressed relative to the contributions from pion loops. 

I now turn to the calculation of finite volume dependence. The loop integrals that appear 
are I\ and I2, proportional respectively to J k G and J k G 2 . They are related by 

J 2 ( m ) = -m^-A-i^m) , (98) 



13 The reader consulting Ref. might be misled by the fact that the leading term in Liischer's formula 
falls like exp(— y / 3/4mX), i.e. more rapidly than exp(— mL). This term, however, is absent for pions because 
it is proportional to the three pion coupling, which is zero. The dominant term falls as exp(— mL), and the 
leading corrections are proportional to exp(— \[2mU). 
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so it is sufficient to consider I\. It is simplest to work in position space 

h{m) = ^ / G(k,m) = ^G A (x = 0,m) . (99) 

/ Ak\<A f 

The subscript A indicates that the position space propagator is smeared over a range 1/A, 
due to the cut-off on the momentum integral. At finite volume we simply substitute Gl (Eq. 
for G. We are interested in the difference 



h(m,L)=I 1 (m,L)-I 1 (m), (100) 

which is given by 

f 2 h(m } L) = Yl G K {n^m) (101) 

« 2 G(n,L M ,m) (102) 

The smearing of the propagator has been dropped in the second line; this is a good approx- 
imation if L >> 1/A, as is true in practice. The precise form of the cut-off is irrelevant. 
Using Eq. (|98|) , the corresponding result for I 2 is 

f 2 T 2 (m,L) = J2 -m 2 ^G(n,L^m) , (103) 

n^O 

In the notation of Gasser and Leutwyler these results are || 

f 2 h(m, L) = 9l (m 2 , Lfj) ; f 2 T 2 (m, L) = m 2 g (m 2 , L,) . (104) 

The finite volume dependence of f n , fx, Bk, By and Ba can now be predicted by substituting 
the results for I\ and I 2 into the formulae presented in previous Sections. 

To actually calculate the /j, it is convenient to use the heat-kernel representation 

G( x ,m) = —— daa- 2 e~ x lAa e~ am . (105) 
(An) 2 Jo 

Using this, the functions go and g\ turn out to be particular cases of the general function 

g r (m,L.) = — — r 'daa r ~ 3 V e -(^W) 2 /^ e -<*™ 2 . (1 6) 
(Air) 2 Jo „ j n 



n^O 



the properties of which have been discussed extensively in Ref. p8 |. The crucial point for 



present purposes is that the sum over n M converges rapidly, as long as mL M >> 1. This is 
because each term falls like exp(— mn^L^), as can be seen by evaluating the a integrals in 
saddle point approximation 



G(x,m) = ^-Jj^Z e — ( 1 + J- + 0(-^) ) (107) 
(Air) 2 y (mx) 6 \ Smx (mx) 2 J 

m 2 -^G(x,m) = -^J— e~ mx \l- — + 0(— ^)) . (108) 
dm 2 (Air) 2 y mx \ 8mx (mx) 2 J 
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To illustrate the properties of the integrals I set L\ = L 2 = L% = L and assume L 4 >> L, 
which is the geometry used in simulations. This means that in the sum over n M only the rij are 
non-zero. It is useful to normalize the results by dividing by the saddle point approximations 
to the terms involving propagation only to the six adjacent periodic images, i.e. rij = (1, 0, 0), 
(-1, 0, 0), (0, 1, 0) etc. These are 



m 2 / ntt 



= 6^^e- (109) 
^n,L) = ^ge- (110) 



Figure |10| shows the ratios I/S, computed numerically, for different numbers of terms in the 



sum over rij. The notation is that Ii(N) includes the first N distances, where "distance" 
means d = y/riiiii. Thus for N = 1 only the six d = 1 terms are kept, while for N = 4, the 
twelve d = v^2, eight d = \/3, arid the six d = 2 terms are also included. Various points are 
noteworthy 

1. It is straightforward to show that the ratios depend only on the product mL. 

2. The results for N — 1 show that the saddle point approximation is accurate to better 



than 10% for mL > 3. If one includes the 1/mL corrections from Eqs. (|107|) and (|108 
then the agreement is better than 1%. 

The full result (N = oo) differs substantially from the leading order term (N = 1) at 
small mL. Typical simulations work down to m n L = 3 — 4, and for these it is clear 



that the leading term is not sufficient. This was not appreciated in Refs. [pT3|, 0, 115 



4. The bulk of the difference is taken up by the N = 2 terms, with N = A being a good 
approximation to the full result. 

5. The ratios for I x and I 2 at N = oo are very similar. 

The need to include not only the d = 1 term in the sum over rij raises the question of 
whether there are other corrections of comparable size. Higher loop diagrams give rise to 
terms which are also proportional to exp(— m n L), but are suppressed by an additional power 
of ml/lQn 2 f 2 . These may be substantial for present lattice pion masses. There are also 
corrections from two loop diagrams suppressed both by a power of m^., and by an additional 
power of exp(— m^L). These are presumably small. 

I close this section with a comment on the effect of the non-zero lattice spacing, a. This 
cuts off the momenta at p pa 2n/a, which should be much larger than the chiral cut-off A. 
This condition is satisfied in present simulations, for which 2n/a > 12 GeV. Such a large 
cut-off should have little effect on the propagator over distances of L or greater. To check 
this, I replace the continuum propagator by that for a lattice scalar field for which the second 
derivative is discretized using the standard nearest neighbor form. This is for illustrative 
purposes only: it is not clear how the pion field should be discretized, or indeed whether 
it should be discretized at all. My choice is really the worst case one can imagine, since 
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discretizations involving more neighbors will give a better approximation to the continuum 
derivative. With the nearest neighbor form, the lattice propagator is 

roc „ 

G(x,y,z,t)= daI x (2a)I y (2a)I z (2a)I t (2a)e- 8a e- am . (Ill) 
J o 

here x, y, z and t are integers giving the distances in lattice units, and I x is the x'th modified 
Bessel function (and should not be confused with the integrals ii and I2 under study). The 
ratio of the integrals using this lattice propagator to those using the continuum propagator 
are shown in Fig. [0]. The results are for N = 1, but there is no significant dependence on 
N . The ratio is close to unity, the difference not exceeding 10% for the values of m n and 
L used in the simulations discussed in the following section. The effect is probably smaller 
than that of higher order chiral corrections. 



7 NUMERICAL PREDICTIONS AND RESULTS 

In this section I compare the predictions for chiral logs and the associated volume dependence 
with numerical results. For the most part the predicted volume dependence is smaller than 
the statistical errors of present simulations. My aim in such cases is to show how close we 
are to testing the predictions. Only for By and Ba is there clear numerical evidence of finite 
volume dependence. 



I use the numerical results from quenched simulations at /3 = 6 presented in Refs. |T3| 
and El. These calculations use staggered fermions, and work on 16 3 x 40 and 24 3 x 40 



lattices. Finite size effects can be studied by comparing results from the two lattices. Using 
the lattice spacing obtained from f w , 1/a = 1.7 GeV, the smaller lattice has L = 1.9 fm, the 
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Figure 11: Ratios of the Ii(N = 1) evaluated using the lattice propagator to those evaluated using the 
continuum propagator. 



larger L = 2.8 fm. Both lattices are probably large enough so that the squeezing of the pion 
wavefunction is a small effect. The lattice pions have masses ranging from m n a = 0.24 to 
0.41, or 0.41 — 0.70 GeV in physical units. The convergence of chiral perturbation theory 
is controlled by the ratio m^/16n 2 f 2 , and this varies from 0.12 to 0.36. Thus higher order 
corrections are likely to be sizable for the heavier pions. 

The pions in these simulations are composed of quarks of masses m q a = 0.01, 0.02 and 
0.03. These are combined in all possible ways, so there are results for pions with both 
degenerate and non-degenerate quarks. The predictions for quenched chiral logs given above 
are only for degenerate quarks. It turns out, however, that the results for f w , B K , By 
and other quantities depend mainly on the average quark mass. When I present results I 
will simply treat all pions as though they contained two degenerate quarks each with the 
average mass. 

Typical chiral logarithmic corrections are numbers of 0(1) multiplied by the integrals I± 
and I<i- The finite volume dependent parts of these integrals are shown in Fig. [12|. These are 
exact results, and are thus a sum of exponentials, multiplied by powers of m n . Nevertheless, 
as the figure shows, the results are well approximated by single exponentials for the range 
of quark masses considered. The magnitude of the li should be compared to the typical 
statistical errors in present numerical simulations, which are rarely as small as 1%. For the 
larger lattice, and for m^a = 0.24 — 0.41, I\ and I2 are always smaller than 0.4%. Thus, for 
practical purposes, L = 24 corresponds to infinite volume. For the smaller lattice, however, 
I2 reaches 3% at m^a = 0.24, which is a level that might be measurable. 

To give an idea of the typical size of finite volume corrections I first consider in a 
theory with four degenerate dynamical quarks. This is the theory represented by a single 
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Figure 12: h and J 2 , for N = oo, using l/o = 1.7 GeV. 



species of staggered fermions in the continuum limit. From Eq. (|24"D we see that 

f n (L)^ f n (L = oo)(l-2h(m n ,L)) . (112) 

It is clear from Fig. [L2] that to observe this effect will require calculations with a precision 
of less than 1%. In the quenched approximation, the prediction is that there is no finite 
size dependence in /„. at leading order. Such dependence enters only due to the diagrams of 
Figs. |^(b) and (c), and is suppressed by m^/167r 2 /^. In fact, the numerical results of Ref. 
[^1 show no evidence for finite size dependence. The errors are large, 5 — 10%, so that this 
does not constitute a test of the predicted difference between full and quenched chiral logs. 

I note in passing that it will be even more difficult to see the volume dependence in 
calculated with dynamical quarks. The result for N degenerate flavors is || 

ml(L)^ml(L = oo)(l + ±T 1 (m 7T ,L)) . (113) 

For the condensate, on the other hand, the situation is slightly better than for f n |J 

(W>)(L) « ®4>KL=oo) (l - ^T x {m^L)) . (114) 

Since the condensate can be calculated more accurately than f n and m n , this equation may 
offer the best chance of testing the predictions of chiral loops in theories with dynamical 
quarks. The only problem is that there is a quadratically divergent perturbative contribution 
to the condensate, to which the prediction does not apply, and which must be subtracted. 

As discussed in Section || quenched chiral logs in m 2 come only from rj' loops. Using 
Eqs. (|89"D and (jS6f), the volume dependence is predicted to be 

ml(L)=ml(L = oo)exp(^^T 2 (m n} L)) , (115) 
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where I have dropped terms suppressed by powers of m\. For N = 3, the phenomenological 
estimate from the physical rf mass gives Ai71q/3 ~ (0.5GeV) 2 . Thus the factor multiplying I2 
ranges from 1.6 to 0.55 for m n = 0.41—0.70 GeV. This means that, although the finite volume 
correction is parametrically enhanced by and is also increased by the exponentiation, 

the actual numerical enhancement is small. The volume dependence is somewhat larger than 
that for the full theory, because it is the integral I2 which appears, rather than I\. The result 
is that, at m n a = 0.24, the correction increases (not its square) by 2.4% for L = 16, 
and by 0.2% for L = 24. The correction decreases rapidly with increasing pion mass. The 
numerical results for show no significant finite size dependence, but the errors are about 



2%, i.e. at same level of the expected effect [p4| . Future quenched calculations may be able 
to test this prediction. 

The form of the volume dependence of the condensate in the quenched approximation 
is the same as that for m 2 , because rn^fl oc m(ipip). Again the numerical results show 
no significant finite volume dependence, but the errors are too large (5-10% for the value 
extrapolated to m q = 0) for this to represent a disagreement with the prediction. 

I now turn to the kaon B-parameters. For degenerate quarks, the predicted volume 



dependence has the same functional form in full and quenched theories. From Eq. (|45[) I 
find 

5B K {L) = B K (L)-B K (L = oo) « B K (L = oo) (h(m K , L) - 2T x (m K , L)) . (116) 
I use irtK rather than to be consistent with the name of the B-parameter; this does not 



imply that the quarks in the "kaon" have different masses. In Fig. [13] I show the fractional 
size of this correction. There is a large cancellation between I2 and I\ in the mass range 
of interest, and the correction does not exceed 0.2%. This is unfortunate as Bk can be 
calculated with errors as small as 1%. The prediction is thus that there be no finite volume 



dependence, and the results of Ref. |L3] are consistent with this. 

As discussed in Section |], the chiral logs in By and B& are enhanced by l/m K over those 
in Bk, and one might hope that the associated finite size effects can be seen. From Eqs. 
(|79D and ©, I find 

B v = \B K - |(/3i + ft)4/2k) - |0*i + £2) - f (71 + l2)h{m K ) + 0(m K ) , (117) 

m K 



where Bk is one-loop corrected result, Eq. (|45|). The enhanced chiral log is that proportional 
to I 2 /m K . The important features of this prediction are: 

• Since l2/m K oc ln(A/mi<-), By should diverge logarithmically in the chiral limit. This 
is unlike Bk, for which the chiral logs, being proportional to m K ln(mK), vanish in the 
chiral limit. 

• The coefficient of this divergence is not related to Bk', it is a new constant proportional 
to /?i + /?2- In vacuum saturation approximation (3i + fa = 4/3. Typically vacuum 
saturation gives the correct sign and order of magnitude; if so, By is predicted to 
diverge negatively. 

• The extent of the enhancement is much greater than that in the correction to the 
quenched m 2 , Eq. (|115|) . The constant /i = m K /(m s + m d ) is about 5 GeV on the 
lattice, roughly ten times the corresponding quantity J Attiq/N. 
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Figure 13: Predictions for the fractional correction SBk(L)/ Bk{L = oo). 



• Since I2 increases with L, By is predicted to become more negative as L decreases. 

Since Ba = Bk — By-, the enhanced log in Ba differs only in sign from that in By. Thus 
Ba is predicted to have a positive logarithmic divergence, and to increase as L decreases. 
Since the numerical results for Bx behave as expected, i.e. they have a smooth chiral limit 
and show no volume dependence, and since the errors in Bk are much smaller than those 
in By and Ba, there is no additional information to be gained from studying Ba- Thus I 
concentrate on By in the following. 



The numerical results of Ref. [IB for By are shown in Fig. [H[ There are two striking 
features: 

• The results appear to be diverging as the chiral limit is approached, with the sign 
expected from the chiral logarithm. 

• There is a significant finite volume dependence, which also has the expected sign. 

It was, in fact, these features that sparked the present study. 
To provide a quantitative test of Eq. ( [L17|) I fit the data using 



By(L) 



16n 2 f 2 x 
H 2 I 2 (m K ,L) 



ci ln(l + l/x) + C2-— h c 3 + c 4 



x 



1 + x 



1 + x 2 



-cr 



:ii8i 



K 



where x = m K /A 2 , A being the cut-off on chiral loop integrals. This function contains the 
chiral logarithm and its associated finite volume dependence, together with a power series 



in m 2 K . The power series is supposed to represent all terms in Eq. (|117|) aside from the 
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Figure 14: Results for By versus m£ in lattice units. The fits are described in the text. 



enhanced chiral logarithm. This includes the normal chiral logs proportional to m 2 K \n{rriK), 
which can be well approximated by a polynomial for the range of masses under consideration. 
I choose the particular form of the first two terms in Eq. (|118|) because the loop integral is 



16vr z 



KA 



G(k, m) 2 = m(l + 1/x) - 1/(1 + x) 



(119) 



The denominator of the C4 term is chosen so that the function is smooth for x ~ 1. Finally, the 
flavor symmetry breaking for staggered fermions is accounted for by applying the prescription 
(|83"D, with Rng 



1, to all terms in Eq. (|118 



Figure [14] shows the best fit. I have taken A to be 0.8 GeV; the precise value is unimpor- 
tant since, to good accuracy, changes can absorbed into the other constants. These constants 
are 

ci = 0.29, c 2 = -1.25, c 3 = 0.63, c 4 = -0.37. (120) 

The fit is reasonable, with a \ 2 °f 1-4 per degree of freedom. I therefore conclude that the 
numerical results are consistent with the expected chiral logarithm. The coefficient c\ is 
smaller than the result expected from vacuum saturation approximation, c\ — 1, but I do 
not consider this to be serious problem for the following reasons. 

• Vacuum saturation can be a poor approximation. For example, the matrix M.k may be 
as much as 50% smaller than the prediction of vacuum saturation ||29|| . What matters 
here, however, are the matrix elements of a pseudoscalar-pseudoscalar operator (cf. 
Eq. (|71~D), and these are likely to be better approximated by vacuum saturation [^0|. 



I have assumed (by setting Rng — 1) that the vertices of the non-Goldstone and 
Goldstone pions are the same, whereas in fact the former could be smaller. The matrix 
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Figure 15: Finite size dependence of By- The line shows the result of the same fit as in Fig. 14. 

elements of some non-Goldstone bilinears are smaller by as much as 25% (Ref. |]20[|). 
Since the loop diagrams involve the matrix element of operators consisting of two 
bilinears, and also a four pion scattering amplitude, a more substantial reduction is 
possible for c\. 

• There could be substantial higher order chiral perturbation theory corrections to the 
coefficient of the enhanced logarithm. That higher order terms can be substantial is 
illustrated by the fact that f^/f 2 = 1.6 for the heaviest lattice kaon. 

The fit does not, however, demonstrate the existence of the chiral logarithm. In particular, 
the coefficient c\ is determined not by the nature of the divergence, but by the finite size 
dependence. If one allows the coefficient of the logarithm to differ from that of the I2 term, 
the fit favors much smaller values for the former. In fact, a fit as good as that displayed 
above is obtained with no logarithmic term at all. This flexibility is possible because x is 
not small, and varies over a small range, 0.5 — 0.8, so the non-logarithmic terms in the fit 
function can mock up a logarithm. 

It is thus the finite volume dependence which is the strongest evidence for quenched chiral 
logarithms. The functional dependence of £?y(16) — .By (24) on tuk is predicted, leaving only 
the overall factor to be determined by the fit. As shown in Fig. [T3|, the fit does a reasonable 
job of representing the dependence on m#. The shape can be altered by higher order chiral 
terms, which are proportional to m 2 K) so one would not expect a perfect match. To provide 
a more rigorous test of the prediction results from smaller kaon masses are needed. 
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8 CONCLUSIONS 



Chiral logarithms can provide powerful consistency checks on numerical simulations. The 
theoretical calculations are well-founded if dynamical quarks are included. I have argued 
that the methods of chiral perturbation theory, suitably modified, can also be applied in the 
quenched approximation. One can imagine a program in which one first tests the results of 
quenched chiral perturbation theory using simple quantities such as /„-, m n and (ipip), and 
then applies the predictions to other quantities in order to test that the physics of pion loops 
is being properly included. Furthermore, as I argued in the Introduction, if the calculations 
are reliable one can use them to estimate the size of the errors introduced by the quenched 
approximation. 

Such a program requires numerical results with errors of 1% or less, calculated on a 
number of lattice sizes at a variety of pion masses. This should be possible in the fairly 
near future, as computer power approaches the TeraFlop milestone. The numerical results 
presented here are encouraging, but not conclusive. What is particularly needed is for the 
calculations to be pushed to smaller pion masses. On a given lattice there is a minimal mass, 
m 7r (min), below which chiral symmetry is restored ||. This minimal mass scales with size 
as m^min) 2 oc L~ 3 . This means that m jr (min)L oc 1/yL, so that, by working on larger 
lattices, one can push to smaller m n L and thereby increase the finite size effects. 

A program of testing finite volume dependence also requires that quenched chiral loga- 
rithms be calculated for other quantities. Examples for which there is hope of good numerical 
results in the next few years include the matrix elements of the AS = 1 operators responsible 
for kaon decays, both CP conserving and violating, and the properties of heavy-light systems, 
such as B and D mesons. Calculations of the finite size effects in fg and fr> in full QCD 



have already been done pi] . For some quantities the quenched calculations appear to be a 
straightforward extension of those presented here. For others, one must face the problem 
of how best to include r( loops. It may be that the perturbative approach of Bernard and 
Golterman [] 1] will be adequate. Or it may be that a resummation, such as that presented 
here the pion mass and the condensate, is required. 

It does not appear to be simple to apply the methods presented here to all quantities. 
In particular, in the calculation of chiral corrections to baryon properties, it may not be 
possible to uniquely assign quark diagrams to chiral perturbation theory diagrams. This 
problem is under study.0 
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